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We consider the superconducting state of d + s symme- 
try with finite concentration of Anderson impurities in 
the limit A s /Ad <c 1. The model consists of a BCS- 
like term in the Hamiltonian and the Anderson impurity 
treated in the self-consistent large- mean field approx- 
imation. 



Increasing impurity concentration or lowering the ratio 
A s /Ad drives the system through a transition from a 
state with two sharp peaks at low energies and exponen- 
tially small density of states at the Fermi level to one 
with N(0) ~ (As/Ad) 2 - This transition is discontinu- 
ous if the energy of the impurity resonance is the small- 
est energy scale in the problem. 
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1 Introduction The order parameter symmetry of 
high-T c compounds is believed to be predominantly of d- 
wave type. However a subdominant s-wave component of 
the order parameter may also occur, especially in materi- 
als with orthorombic distortions. One of such compounds 
is YBCO, exhibiting strong structural distortion and a 
substantial anisotropy in the London penetration depth in 
the a-b plane. (TJ Raman scattering J2] provided evidence 
for a 5% admixture of s-wave component while thermal 
conductivity measurements in rotating magnetic field [3] 
placed an upper limit of 10%. Angle-resolved photoemis- 
sion spectroscopy JU on monocrystalline YBCO gave the 
ratio of 1.5 for gap amplitudes in the a and b directions in 
the CuOi- Measurements of a Josephson current between 
monocrystalline YBa2Cu307 and s-wave Nb showed that 
the obtained anisotropy could be explained by a 83% 
(i-wave with a 17% s-wave component. Q In another ex- 
periment on YBCO/Nb junction rings the s to d gap ratio 
in optimally doped YBCO was estimated to be 0.1.151 

Inelastic neutron scattering on monocrystalline and un- 
twinned samples of YBCO lead to magnetic susceptibili- 
ties with intensities and line shapes breaking the tetragonal 
symmetry.! 7,8 ,9. 10| It was shown that these data may be 
interpreted within an anisotropic band model with an order 
parameter of mixed d and s symmetry. ifTTl 

Superconducting states with mixed symmetry are also 
considered in other classes of compounds, e.g. in the re- 
cently discovered ferropnictides. lfl2l 



The effects of dilute concentrations of magnetic and 
nonmagnetic point defects on a BCS superconductor of 
pure s or rf-wave symmetry were intensively studied in the 
past and are well known. In a ci-wave superconductor with 
lines of order parameter nodes any amount of disorder in- 
duces a nonzero density of states at the Fermi level. In an s- 
wave system only magnetic impurities change the response 
of the superconducting state. For sufficiently strong cou- 
pling between the impurity and the conduction band bound 
states may appear in the energy gap. lTOl 

In an earlier work on nonmagnetic impurities inad + 
s-wave superconductor it was shown that in the unitary 
limit a nonzero density of states (DOS) at the Fermi level 
appears above certain critical impurity concentration, de- 
pending on the size of an s-wave component. lfT4l However 
the low-energy DOS is mostly featureless since the s-wave 
component prevents a buildup of states due to nonmagnetic 
scattering. In contrast the presence of magnetic impurities 
in a superconductor with an s-wave component may result 
in sharp peaks in the low-energy DOS for small concentra- 
tion of defects, provided the energy scale associated with 
the impurity resonance is small. 

2 Model and Results We consider the order param- 
eter of d + s symmetry on a cylindrical Fermi surface. 
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where A s and Ad are amplitudes of s- and ri-wave compo- 
nent respectively. We assume 9 = and A s <C Ad- 

The superconductor is treated in a BCS approximation. 
The magnetic scatterer is modelled as an Anderson impu- 
rity treated within the slave boson mean field approach. 1 1 5 1 
The low energy physics is dominated by the presence of 
strongly scattering impurity resonance. The self-consistent 
self-energy equations describing the interplay between the 
superconducting and magnetic degrees of freedom have the 
following form, 



nN 

LU = LO + - — —r 
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where u5(u>), A(A) is the renormalized frequency and or- 
der parameter of conduction electrons (impurity) respec- 
tively. 

We assume Ad/ D = 0.01, where 2D is the bandwidth 
of the conduction electron band. In the equations above r 
is the hybridization energy between the impurity and the 
conduction band and e/ is the resonant level energy. Here 
we assume constant density of states in the normal state 
No = 1/2D and do the calculations for a nondegenerate 
impurity, N = 2. Brackets denote average over the Fermi 
surface. 

Initial results for the density of states of this system 
were presented in an earlier paper. lfT6l For small impurity 
concentration n, such that A s and Ad are not significantly 
affected, there are two peaks located symmetrically near 
the gap center, provided e/ <C -T <C A s . For larger n 
the two peaks merge into one peak centrally located at the 
Fermi energy. 

Fig. 1 shows the dependence of the density of states 
N(Q) at Ep as a function of impurity concentration. 
At small n, N(0) is exponentially small, N(0)/N ~ 
(A s /Ad) cxp(— a(A s /Ad) 2 /n), where a is a numerical 
factor. The critical concentration no for the discontinuous 
transition to N(Q)/No ~ A s /Ad is a quadratic func- 
tion of A s /Ad- For n > no, N(Q) approaches DOS of a 
d-wave superconductor in both the magnitude and its func- 
tional dependence on n. These relations are valid when 
e f < r,A s . 

Qualitatively similar scaling of N(0) as a function of 
A s /Ad in presence of nonmagnetic impurities in the uni- 
tary limit was obtained in ref. Ifl4ll . The difference between 
magnetic and nonmagnetic defetcs in a superconductor 
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Figure 1 Logarithm of the density of states at the Fermi 
level as a function of impurity concentration for several 
values of the ratio A s j Ado- The resonant impurity level e/ 
is close to the Fermi level. r/D is fixed at 0.001. 
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Figure 2 Position of the peak in the conduction electron 
DOS as a function of impurity concentration for d + s su- 
perconductor for several values of the resonant level en- 
ergy. Energy is scaled by half of the conduction electron 
bandwidth D and r/D = 0.001. The suppression of A s 
and Ad was not taken into account. 



with A s / 'Ad <C 1 is the character of the transition from 
the exponentially small A^(0) to finite A^O) in the strong 
scattering limit. In contrast to nonmagnetic impurities, the 
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transition caused by resonantly scattering magnetic impu- 
rities is discontinuous, and there are two sharp peaks of 
N(co) on both sides of the transition. lfT6l The position of 
peaks as a function of impurity concentration for differ- 
ent values of is shown in Fig 2. These peaks strongly 
alter the low-energy and low-frequency response and may 
be detected in thermodynamic or transport measurements. 
The increase of A s shifts the resonances towards Ep and 
makes them more narrow. 

The d-wave component of the superconducting order 
parameter is more sensitive to pair breaking by magnetic 

impurities in the limit T <C Tk, where T K = J ' T 2 + e 2 , 

than the s-wave part. Depending on the relative size of Ad 
and A s there may be another impurity transition for larger 
n, when the order parameter nodes disappear due to van- 
ishing of the d-wave component and a full gap opens up. 
The impurity peak is then split and (0) falls to zero again. 
Hower a detailed description of this possibility requires a 
careful analysis involving four energy scales: Ad, A s , T, 
and e / . 

Can such transition be observed experimentally? The 
transition may be tuned either by varying impurity con- 
centration or the ratio A s j Ad- The density of states of the 
normal state at Ep, No, also has an effect. It appears in 
equations|2]and[3] One should bear in mind, however, that 
both the superconducting transition temperature and the 
impurity resonance energy scale Tk are exponentially sen- 
sitive to changes of No. Experiments conducted at very low 
temperatures may give different signatures of low-energy 
behavior depending on the location of the system on the 
phase diagram relative to the impurity critical point. 

In a d+s-superconductorin the limit of high concentra- 
tion of point defects the d-wave component vanishes while 
the s-wave part of the order parameter remains unaffected. 
The situation is different in presence of magnetic scatter- 
ers. If the s-wave component is small and the energy scale 
of the resonance due to impurity scattering is at most of 
the order A s , increasing impurity concentration may drive 
A s — > while Ad will be reduced and finite. This follows 
from the fact that the largest pair breaking occurs when 
the energy scales of the impurity resonance and the order 
parameter are comparable. If Ad 3> Tk, the rate of sup- 
pression of Ad with increasing n is small. 

The d + s superconductor with a subdominant s-wave 
component doped with magnetic impurities has a differ- 
ent phase diagram in the T-n plane compared to the same 
superconductor with nonmagnetic defects. If the impurity 
resonance scale Tk is comparable to A s and A s <C Ad, 
the s-wave component may vanish at T ~ Tk- 

3 Conclusions The T = impurity transition dis- 
cussed in this work may be detected in the low temperature 
limit. While the slave boson mean field formalism used in 
this paper cannot be applied at temperatures exceeding Tk, 
we may also qualitatively describe the expected behavior 
of the system as a function of temperature. At finite but 



small impurity concentration there may be even four phase 
transitions as a function of temperature: normal to d-wave, 
d-wave to d + s-wave, d + s to d-wave, and d-wave to 
d + s-wave. Due to the difference of magnitudes the small 
A a may be driven to zero faster with increasing n than 
Ad- Nonmonotonic behavior around T ~ Tk is a conse- 
quence of strong scattering by the resonance state forming 
on the impurity site. As T — > 0, the impurity scattering 
becomes weaker and the s-wave component may appear 
again. Whether this particular scenario is realized, depends 
on the relative size of energy scales: A s /Ad, A s /Tk, and 

tf/r. 
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